Many-body singlets by dynamic spin polarization 
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We show that dynamic spin polarization by collective raising and lowering operators can drive a spin ensemble 
from arbitrary initial state to many-body singlets, the zero-collective- spin states with large scale entanglement. 
For an ensemble of N arbitrary spins, both the variance of the collective spin and the number of unentangled 
spins can be reduced to O(l) (versus the typical value of O(N)), and many-body singlets can be occupied 
with a population of ~ 20% independent of the ensemble size. We implement this approach in a mesoscopic 
ensemble of nuclear spins through dynamic nuclear spin polarization by an electron. The result is of two-fold 
significance for spin quantum technology: (1) a resource of entanglement for nuclear spin based quantum infor- 
mation processing; (2) a cleaner surrounding and less quantum noise for the electron spin as the environmental 
spin moments are effectively annihilated. 

PACS numbers: 76.70.Fz,42.50.Dv,03.67.Bg,71.70.Jp 



Many-body singlets (MBS) are the zero-collective- spin 
states of a spin ensemble with large scale quantum entangle- 
ment and zero spin uncertainties. They appear in a variety of 
contexts in quantum physics and in condensed matter physics, 
e.g. as horizon states of the quantum black hole ll|], and as 
ground states of quantum antiferromagnetic models [2] . Their 
special characteristics place them at the center of attention for 
quantum applications. First, MBS are invariant under a simul- 
taneous unitary rotation on all spins. This makes MBS suit- 
able for spanning a decoherence-free subspace [3], for quan- 
tum communications without a shared reference frame 0], 
and for metrology of the spatial gradient or fluctuations of ex- 
ternal fields Q]. Second, MBS is an extreme example for the 
squeeze of spin uncertainties The collective spin has 

zero variance in all directions and thus a source of quantum 
noise is removed, e.g. in the context of a quantum object af- 
fected by a spin bath. Third, MBS contain large scale quantum 
entanglement: every spin is entangled with the rest part of the 
ensemble. An example of a pure MBS is the product of two- 
qubit singlets (Bell pairs). In the maximally mixed state of all 
MBS, the distillable bipartite entanglement is logarithmic in 
the ensemble size . 

Despite the successful generation of photonic analog of 4- 
qubit singlets by parametric down conversion flfioll. realiza- 
tion of MBS in a general spin ensemble is an outstanding goal 
awaiting technically feasible approaches. Theoretical study 
shows that spin squeezing based on quantum non-demolition 
measurement can reduce the total collective spin variance of 
an atomic ensemble by a factor of 5 in the lossless case Q]. 
However, in such squeezed state the weighting of MBS is 
small and vanishes in large TV limit. 

Here we introduce a conceptually new approach for squeez- 
ing of collective spin uncertainties and generation of large 
scale entanglement. The approach uses collective spin raising 
and lowering operations only, and is applicable to an ensem- 
ble of N arbitrary spins initially on arbitrary state. The state 
after squeezing is significant in figures of merit: in the low 
loss limit, MBS are occupied with an TV-independent popula- 
tion of ~ 20%, and both the variance of the total collective 



spin and the number of spins unentangled with the rest are 
0(1) (versus the typical values of O(N)). We implement this 
approach in a mesoscopic ensemble of nuclear spins, a spin 
system of extensive interests either as a noise source or as a 
superior information storage in quantum technology. The im- 
plementation uses only generic features of dynamic nuclear 
spin polarization processes by an electron, and is applicable 
to various electron-nuclear spin systems. Distillation of MBS 
can be realized by post-selection based on measurement of 
the electron spin. MBS can be a valuable resource of quan- 
tum entanglement for nuclear spin quantum information pro- 



cessing HI lMl3h . In electron spin based quantum computation 



schemes, preparing the peripheral nuclear spins into MBS re- 
sults in a cleaner surrounding and hence improved quantum 
coherence of the electron spin. 

We refer to the definition of spin squeezing in the general- 
ized sense Q,@], where the degree of squeezing is quantified 
by (J 2 ), with J = J2n=i In being the total collective spin 
for an ensemble of N particles with equal or different spins. 
(J 2 ) = indicates perfect squeezing where the TV spins are 
in the MBS. (J 2 ) gives an upper bound on the number 

of spins unentangled with others where s is the average spin 
per particle |a,|2D. States of the spin ensemble can be grouped 
into multiplets, i.e. irreducible invariant subspaces of the to- 
tal spin. A multiplet { | J, M, a 1 }) , M = - J, . . . , J} will be 
denoted in short as { J, a 1 }}, where a 1 } is a general index for 
distinguishing the set of orthogonal (2 J+l)-dimensional mul- 
tiplets. The aim is to transfer population from all multiplets to 
those singlets with J = 0. 

Key to our squeezing approach is to apply raising operator 
of the form j\ — jg on the spin coherent states \J,—J,ctj). 
Here the ensemble is partitioned arbitrarily into two subsets A 
and B with collective spin ja and js respectively (Fig. 1(a)). 
We find the key identity 

(J + l,-J+l,a J+ i|(j'+-j'+)|J,- 

(J,-J,aj\ 0a~ Jb) \ j + 

= -[(J+l)(2J+l)]-*. 

Moreover, (J', M', aj, \ j+ - j+ | J, 
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FIG. 1: (a) The red dashed line partitions the spin ensemble into 
two subsets A and B with collective spin ja and jb respectively. 
The green dashed line gives a different partition of the ensemble into 
subset C and D. (b) An example of how the operator j\ — cou- 
ples various basis states | J, M,2a,2b) (solid arrows). The absolute 
value squared of the transition matrix elements are indicated with 
the thickness of the arrows. The transitions related to the spin co- 
herent states | J, — J) are highlighted. Hollow vertical arrows show 
the coupling by J~ . (c) Schematics of the population transfer rates 
between multiplets under the condition that each multiplet is initial- 
ized on the spin coherent state | J, — J) when the j\ — jg operator 
is applied. The transfer rate of { J, aj} — > {J,aj } is identical to 
the rate of the backward transfer {J, a 1 } } — > {J, a 1 }}. The rate of 
{J + 1, a q J+1 } {J, a 1 }} is by a factor of (J + 1)(2 J + 1) larger 
than that of the backward transfer {J, a 1 }} —> {J + 1, 



J | > 1 or M' ^ — J + 1. Thus, under the condition that each 
multiplet is initialized on the spin coherent state, application 
of the j\ —jg operator tends to transfer populations from mul- 
tiplets of larger dimension to multiplets of smaller dimension 
(Fig. 1(c)). The transfer rate of {J + 1, -> {J, «j} is 

by a factor of(J+l)(2J+l) larger than that of the backward 
transfer {J, a 1 }} { J + 1, a q J+1 }. 

Squeezing of collective spin uncertainties can therefore 
be realized by dynamic spin polarization with the lower- 
ing operator J~ and raising operators of the form j\ — 
jg. Consider the use of two such operators j\ — j% 
and 2c ~ 3d wnere C an d D constitute a different bi- 
partition of the ensemble (Fig. 1(a)). The Hilbert space 
can be divided into independent subspaces according to 
the quantum numbers {jAnD, 3bhd, Jbhc, 3Anc} conserved 
by the raising/lowering operations, and their values deter- 
mine the number of (2J+l)-dimensional multiplets n(J). 
If J~ is applied more frequently such that the system 
is in spin coherent states every time j\ ~ 3 b or 3c ~ 




FIG. 2: (a) An example of the squeezing control by dynamic spin 
polarization with the operators J - , j\ ~ 3b m ^ 3c ~ 3d- 5 ) 
and (c) show simulation results of the control in (a). The sys- 
tem is initially in the completely mixed state in the subspace with 
{Jahd = 7/2JBHD = 7/2J B nc = 7/2J A nc = 7/2}. The 
polarization rates: 10 _3 A^ = A G , and the delay time r = 2/A . 
(b) p( J, 2 a, 2 b) gives the population of the multiplet {J, ja, 2 b} in 
the state at t = At. (c) P(J) = J2j A ,j B p(JJaJb) gives the in- 
tegrated probability of finding the system with total collective spin J 
at various time. At t = 5r, MBS are occupied with the population 
P(J = 0) = 0.21, and (J 2 ) = 2.42. 



jp is applied, we find the steady-state in each subspace: 

P = Ejf(J)E n k L J ? \J,-J,a k j){J,-J,a k j\ where /(J) = 
(J + 1)(2J + l).f (J + 1). Most subspaces contain at 
least one MBS Q, and n (J) < ra(0)(2J + l). Thus, 
in the steady state, MBS are occupied with a population 
n(0)/(0) > Ej^(J)]" 1 = 0.20, and the variance (J 2 ) < 
[Y,jg(J)Y 1 Y,jHJ + ^g(J) = 2.44, where g(J) = 
(2J+l)[n£o 1 (i + l)(2i + l)] \ 

Dynamic spin polarization by the collective raising and 
lowering operators can be described by Lindblad terms in 
the master equation p = -\ ELi(4^mP + pL^Lm - 
2L m pLl n ) where L x a/A^J", L 2 V^Qa 3b)> 
and L 3 = \/~A~ (jc - 3d)' The ^~ and 3 a ~ 3b operators 
are applied with the rates A^ | (ipf \ J~ \ ) | 2 and A \(i/>f \j\ — 
3b 1^)1 2 respectively, and the squeezing scheme requires 
that the former rate shall always be larger. We note that 

(J,MJ a Jb\0a ~3b)0a ~ 3b)\ j ^ m ^3aJb) increases 
with the decrease of J and reaches the maximal value of 
~ (J a + 3b? for small J, while (J, M\J + J~\J, M) - J 2 . 
Thus we find the requirement A^/A > (jAnD + 3bc\d + 
3Bnc + jAnc) 2 , the latter quantity ~ ANs 2 in an ensemble 
of N spin s particles. Spin decoherence causes the population 
decay of MBS with a rate ~ Nj n with j n being the single 
spin decoherence rate. The low loss condition is therefore de- 
fined as 



4Ns' 



-Ah > A Q ^> 7 n , where spin decoherence has 



negligible effect on the squeezing efficiency 11511 . 



We numerically demonstrate a squeeze control where j\ — 
jg and — j+ are applied in alternating fashion (see 
Fig. 2(a)). With spin decoherence neglected under the low loss 
condition, calculation can be significantly simplified for this 
choice of control and a moderately large spin system can be 
simulated. In the interval when J~ and j\ — jg (or — j+) 
are applied, the relevant Hilbert space can be further divided 
into independent subspaces according to the quantum num- 
bers {jaiJb} (or {jc, 3d})- If the duration r of each interval 
is sufficiently large to ensure the reach of steady state, we sim- 
ply need to solve for steady state in each small subspace and 
keep track of the basis transform between {| J, M, ja, 3 b)} 
and {| J, M, jc, 3d)} upon the switch of raising operators. 
Off-diagonal coherence is found to be negligible which can 
further simplify the calculation. An example of the simulated 
squeezing dynamics is given in Fig. 2(b-c). The initial den- 
sity matrix is the completely mixed one in the subspace with 
{jAnD = 7/2 J BnD = 7/2J BnC = 7/2J A nc = 7/2}. 
The polarization rates are 10 _3 A/ l = A , and r = 2/A G . 
Fig. 2(b) shows that the population distribution p(J, j^, 3b) 
among the multiplets indeed approaches the steady state value 
after a time of 4r. At t = 5r, MBS are occupied with a pop- 
ulation of 0.21 and (J 2 ) = 2.42. 

In an ensemble of nuclear spins, the applications of the two 
types of operators J~ and j\ — are realized in the pro- 
cess of dynamic nuclear spin polarization (DNSP), a major 
tool for manipulation of nuclear spins 1116142911 . We consider 
the hyperfine interaction Hq = ^2 n \^ (r n )\ 2 I n ■ A • S cou- 
pling the electron spin S to peripheral lattice nuclear spins 
I n . A is the hyperfine coupling constant in tensor form, and 
the position dependence of coupling enters through the enve- 
lope function ip (r) of the electron only. Hq describes gen- 
erally the hyperfine interaction of electron or hole system in 
quantum dots or shallow donors formed in group IV or III- 
V materials JIolEll]. In most DNSP schemes, Hq induces 
the electron-nuclear flip-flop in passing electron spin polar- 
ization to the nuclei and the energy cost is compensated by 
emission/absorption of phonons or photons 1119142311 . These 
DNSP schemes are termed as the dc type hereafter. Alterna- 
tively, DNSP can also utilize the ac correction to the hyperfine 

coupling: H ac = V 

cos(ut) I n • A • S 

when an ac electric field induces an electron displacement 
d u cos(cjt), with energy cost for electron-nuclear flip-flop di- 
rectly supplied by ac field I27U29I1 . Such DNSP process is 
termed hereafter as the ac type. 

For nuclear spins on the periphery of an electron, MBS 
can be realized by combining dc and ac DNSP processes 
which polarize nuclear spins in opposite directions with the 
operators Y, n \^ On)| 2 K and E n ^ \^ ( r n)\ 2 In respec- 
tively. Here fi is the direction of ac electric field. The lat- 
tice sites with equal electron density |^(r)| 2 are grouped 
into coordination shells. On each shell, \ip (r n )| 2 I~ and 

En ^7 1^ On) | 2 In are of the character of J~ and j\ - j% 
respectively. Under the influence of incoherent electron spin 
dynamics in the DNSP process, the large shell-to- shell differ- 
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FIG. 3: (a) Schematics of an electron in a nuclear spin bath. The first 
coordination shell has 4 nuclear spins (green color) and the second 
shell has 8 nuclear spins (blue color), (b) Upper part: dc hyperfine 
coupling coefficients at the various lattice sites. Lower part: ac hy- 
perfine coupling coefficients with ac electric field in x direction. The 
heights of the bar give the magnitude, (c) The ac hyperfine coupling 
decomposed into two terms H ac = il\ c + &ac- — &i0a — 
j+)S-e^E x + c.c, and H 2 ac = ~a 2 {ji ~ J^S'e^E, + c.c. 
The 4 sites with positive coupling coefficients in the upper part form 
subset A and the 4 sites with positive coupling in the lower part form 
subset C, and B (D) is the complement of A (C). (d) Schematic 
of the DNSP control where the switching between dc and ac DNSP 
is concatenated with the switching of the ac electric field between x 
and y directions. 



ence in the dc hyperfine coupling strength causes loss of inter- 
shell coherence in a timescale much faster than the squeez- 
ing. Thus different coordination shells can be independently 
squeezed towards MBS. 

Fig. 3(a) shows the schematic of an electron with a 2D 
Gaussian envelope function. The 12 lattice nuclear spins form 
two coordination shells according to the dc hyperfine cou- 
pling strength (Fig. 3(b)). For the green shell with 4 lattice 
nuclear spins, H ac = aS~(j\ - j^)e iuJt E x + aS~(j^ - 
j^)e tuJt E y + c.c, where E x ^ is the ac electric field in the 
x (y) direction. Fig. 3(d) shows the schematic of the DNSP 
control where the switching between dc and ac DNSP is con- 
catenated with the switching of the ac electric field between x 
and y directions. Numerical simulation of such a control has 
been given in Fig. 2. The blue shell represents the more gen- 
eral case where nuclear spins are polarized in the ac DNSP 
process by a\(j\ — 3b) + ^Uc ~ Jd)* a nnear superposi- 
tion of raising operators of the desired form (Fig. 3(c)). The 
identity in Eq. (Q]) obviously holds if j\ — jg is replaced 
by — jf+, thus we have this same identity for their linear 
superposition as well. Numerical simulations confirm that op- 
erators of this new form are equally efficient in the squeezing 
effect iH. 

Interaction between neighboring nuclear spins causes spin 
diffusion and spin dephasing which can result in loss of MBS. 
The dipolar interaction between neighboring lattice sites is of 
the strength ~ 10 Hz. Nuclear spin diffusion by the 
coupling terms is efficiently suppressed when the shell-to- 
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shell inhomogeneity in the hyperfine coupling is large. By the 
In^m term, the nuclear spins are subject to a dipolar magnetic 
field dependent on the configuration of their neighbors, which 
leads to dephasing with a rate j n ~ 10 — 100 Hz. To realize 
efficient squeezing, fast DNSP mechanisms are desired. 

For optically controllable electron spin, e.g. in quantum 
dot or impurity in III-V semiconductors, fast dc DNSP can be 
realized by the hyperfine-mediated optical excitation of spin- 



forbidden excitonic transitions 112 lL 13211 . Assuming the elec- 
tron Zeeman splitting uo e ~ 0.2 GHz, the intrinsic broadening 
of charged exciton j t ~ 0.2 GHz, and an optical Rabi fre- 
quency ft ~ 3 GHz for the excitonic transition, we estimate 
the DNSP rate: = 9l 4^~ ~ 10 MHz on a coordination 
shell with hyperfine coupling a = 3 MHz. For other electron- 
nuclear spin system, fast dc DNSP may be realized through 
the bath-assisted electron-nuclear flip-flop in the presence of 
an efficient energy dissipation channel, e.g. an electron Fermi 
sea in nearby leads I120ll . 

~2 

ac DNSP is of the rate A = — where 7 S is the broad- 
ening of the electron spin resonance 112 811 . The magnitude of 
the ac hyperfine interaction, a, depends on the strength of ac 
electric field and the inhomogeneity of the electron envelop 
function. Giving the phosphorus donor in silicon as an exam- 
ple, we have the first several shells: (A, 6.0, 6), (B, 4.5, 12), 
(C, 3.3, 4), (D, 2.2, 12) and (F, 1.7, 12) where the first letter 
is the label of the shell by convention, the second number is 
the hyperfine coupling strength in unit of MHz, and the third 
is the number of equivalent sites on the shell [33]. The dis- 
tance between neighboring shells is in the order of ~ 0.1 nm. 
Thus, we estimate a ~ MHz by a moderate displacement of 
the electron d u ~ 0.1 nm. Assuming 7 S ~ 0.1 GHz, A can 
be of ~ 10 kHz on these shells. Thus, we conclude that the 
low loss condition can indeed be satisfied for nuclear spins on 
the periphery of a strongly confined electron. 

The work was supported by the Research Grant Council 
of Hong Kong under Grant No. HKU 706309P. The author 
acknowledges helpful discussions with L. J. Sham, H. Y. Yu 
and X. D. Xu. 



[1] E. R. Livine and D. R. Terno, Phys. Rev. A 72, 022307 (2005). 
[2] G. Misguich and C. Lhuillier, Frustrated Spin Systems (World 

Scientific, 2004), chap. 5, p. 229. 
[3] M. Bourennane, M. Eibl, S. Gaertner, C. Kurtsiefer, A. Cabello, 

and H. Weinfurter, Phys. Rev. Lett. 92, 107901 (2004). 
[4] S. D. Bartlett, T. Rudolph, and R. W. Spekkens, Phys. Rev. Lett. 

91, 027901 (2003). 
[5] G. Toth and M. W. Mitchell, New J. Phys. 12, 053007 (2010). 
[6] D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore, and 

D. J. Heinzen, Phys. Rev. A 46, R6797 (1992). 



U 
[8 

[9 

[10 



[11 
[12 



[13 
[14 



[15 
[16 



[17 
[18 
[19 



[20 
[21 
[22 



[23 



[24 
[25 



[26; 

[27 

[28; 
[29 
[30 

[31 

[32; 

[33 



M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

A. S. Sorensen and K. Molmer, Phys. Rev. Lett. 86, 4431 
(2001). 

G. Toth, C. Knapp, O. Guhne, and H. J. Briegel, Phys. Rev. 
Lett. 99, 250405 (2007). 

M. Eibl, S. Gaertner, M. Bourennane, C. Kurtsiefer, 
M. Zukowski, and H. Weinfurter, Phys. Rev. Lett. 90, 200403 
(2003). 

B. E. Kane, Nature 393, 133 (1998). 

M. V. Gurudev Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, 
F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, Sci- 
ence 316, 1312 (2007). 

P. Neumann, N. Mizuochi, F. Rempp, P. Hemmer, H. Watan- 
abe, S. Yamasaki, V. Jacques, T. Gaebel, F. Jelezko, and 
J. Wrachtrup, Science 320, 1326 (2008). 
The states unconnected with MBS are of ~ 1% of the entire 
Hilbert space if two raising operators of the form j\ — jg are 
used, and can be reduced to ~ 0.01% if three such operators 
are used in the dynamic spin polarization. 

H. Y. Yu and W Yao, unpublished. 

A. Greilich, A. Shabaev, D. R. Yakovlev, A. L. Efros, I. A. Yu- 
gova, D. Reuter, A. D. Wieck, and M. Bayer, Science 317, 1896 
(2007). 

D. J. Reilly, J. M. Taylor, J. R. Petta, C. M. Marcus, M. P. Han- 
son, and A. C. Gossard, Science 321, 817 (2008). 

X. Xu, W Yao, B. Sun, D. G. Steel, A. S. Bracker, D. Gammon, 

and L. J. Sham, Nature 459, 1105 (2009). 

A. Tartakovskii, T. Wright, A. Russell, V. Falko, A. Vankov, 

J. Skiba-Szymanska, I. Drouzas, R. Kolodka, M. Skolnick, 

P. Fry, et al., Phys. Rev. Lett. 98, 026806 (2007). 

J. Danon and Y. Nazarov, Phys. Rev. Lett. 100, 056603 (2008). 

V. Korenev, Phys. Rev. Lett. 99, 256405 (2007). 

A. S. Bracker, E. A. Stinaff, D. Gammon, M. E. Ware, J. G. 
Tischler, A. Shabaev, A. L. Efros, D. Park, D. Gershoni, V. L. 
Korenev, et al., Phys. Rev. Lett. 94, 047402 (2005). 

B. Urbaszek, P. Braun, T. Amand, O. Krebs, T. Belhadj, 
A. Lemaitre, P. Voisin, and X. Marie, Phys. Rev. B 76, 201301R 
(2007). 

K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 (2004). 

C. Latta, A. Hogele, Y. Zhao, A. N. Vamivakas, P. Maletinsky, 
M. Kroner, J. Dreiser, I. Carusotto, A. Badolato, D. Schuh, 
et al., Nat. Phys. 5, 758 (2009). 

I. T. Vink, K. C. Nowack, F. H. L. Koppens, J. Danon, Y V. 
Nazarov, and L. M. K. Vandersypen, Nat. Phys. 5, 764 (2009). 

E. Laird, C. Barthel, E. Rashba, C. Marcus, M. Hanson, and 

A. Gossard, Phys. Rev. Lett. 99, 246601 (2007). 

M. Rudner and L. Levitov, Phys. Rev. Lett. 99, 246602 (2007). 
E. Rashba, Phys. Rev. B 78, 195302 (2008). 

B. Eble, C. Testelin, P. Desfonds, F. Bernardot, A. Balocchi, 
T. Amand, A. Miard, A. Lemaitre, X. Marie, and M. Chamarro, 
Phys. Rev. Lett. 102, 146601 (2009). 

J. Fischer, W Coish, D. Bulaev, and D. Loss, Phys. Rev. B 78, 
155329 (2008). 

E. A. Chekhovich, M. N. Makhonin, K. V. Kavokin, A. B. 
Krysa, M. S. Skolnick, and A. I. Tartakovskii, Phys. Rev. Lett. 
104, 066804 (2010). 

E. B. Hale and R. L. Mieher, Phys. Rev. 184, 739 (1969). 



